Exact results for SU(3) spin chains: 
trimer states, valence bond solids, and their parent Hamiltonians 
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We introduce several exact models for SU(3) spin chains: (1) a translationally invariant parent 
Hamiltonian involving four-site interactions for the trimer chain, with a three-fold degenerate ground 
state. We provide numerical evidence that the elementary excitations of this model transform under 
representation 3 of SU(3) if the original spins of the model transform under rep. 3. (2) a family of 
parent Hamiltonians for valence bond solids of SU(3) chains with spin reps. 6, 10, and 8 on each 
lattice site. We argue that of these three models, only the latter two exhibit spinon confinement 
and a Haldane gap in the excitation spectrum. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.10.Dg, 32.80.Pj 
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Beginning with the invention of the Bethe Ansatz in 
1931 [1] as a method to solve the 5 = 5 Heisenberg chain 
with nearest neighbor interactions, a significant share of 
the entire effort in condensed matter physics has been de- 
voted to the study of quantum spin chains. Faddeev and 
Takhtajan [2] discovered in 1981 that the elementary ex- 
citations (now called spinons) of the spin-1/2 Heisenberg 
chain carry spin 1/2 while the Hilbert space is spanned 
by spin flips, which carry spin 1. The fractional quanti- 
zation of spin in spin chains is comparable to the frac- 
tional quantization of charge in quantized Hall liquids [3] . 
In 1982, Haldane [4] identified the 0(3) nonlinear sigma 
model as the effective low-energy field theory of SU(2) 
spin chains, and argued that chains with integer spin pos- 
sess a gap in the excitation spectrum, while a topological 
term renders half-integer spin chains gapless. 

The general methods — the Bethe Ansatz and the use of 
effective field theories including bosonization — are com- 
plemented by a number of exactly solvable models, most 
prominently among them the Majumdar-Ghosh (MG) 
Hamiltonian [5] for the 5 = \ dimer chain, the AKLT 
model [6] as a paradigm of the gapped 5=1 chain, and 
the Haldane-Shastry model (HSM) [7-9]. In the HSM 
the wave functions for the ground state and single-spinon 
excitations are of a simple Jastrow form, elevating the 
conceptual similarity to quantized Hall states to a for- 
mal equivalence. One of the unique features of the HSM 
is that the spinons are free in the sense that they only 
interact through their half- Fermi statistics [10-12]. The 
HSM has been generalized from SU(2) to SU(n) [13-16]. 

For the MG and the AKLT model, only the ground 
states are known exactly. Nonetheless, these models have 
amply contributed to our understanding of many aspects 
of spin chains, each of them through the specific concepts 
captured in its ground state [17-24]. 

In the past, the motivation to study SU(n) spin sys- 
tems with 11 > 2 has been mainly formal. The Bethe 
Ansatz has been generalized to multiple component sys- 
tems by Sutherland [25], and has been applied to the 
SU(n) HSM [13, 14]. The effective field theory descrip- 
tion of Haldane yielding the distinction between gapless 
half-integer spin chains and gapped integer spin chains, 



however, cannot be directly generalized to SU(n) chains, 
as there is no direct equivalent of the CP 1 representation 
used in Haldane's analysis. 

In recent years, ultracold atoms in optical lattices 
have provided a framework for model realizations of var- 
ious problems of condensed matter physics, including 
the phase transition from a superfluid to a Mott insu- 
lator [26, 27], the fermionic Hubbard model [28, 29], and 
SU(2) spin chains [30, 31]. In particular, the Hamiltoni- 
ans for spin lattice models may be engineered with polar 
molecules stored in optical lattices, where the spin is rep- 
resented by a single-valence electron of a heteronuclear 
molecule [32] . Systems of ultracold atoms in optical lat- 
tices may further provide an experimantal realization of 
SU(3) spin systems, and in particular antiferromagnetic 
SU(3) spin chains, in due course. A simple and intrigu- 
ing possibility is to manipulate an atomic system with 
total angular momentum F = | such that it simulates 
an SU(3) spin. For such atoms, one effectively suppresses 
the occupation of one of the "middle" states, say the 
F z = — i state, by shifting it to a higher energy while 
keeping the other states approximately degenerate. At 
sufficiently low temperatures, one is left with three inter- 
nal states F z = — §,+i,+f, which one identifies with 
the colors "blue" , "red" , and "green" of an SU(3) spin. 
In leading order, the number of particles of each color 
is now conserved, as required by the SU(3) symmetry. 
If one places one of these atoms at each site of an opti- 
cal lattice and allows for a weak hopping, one obtains an 
SU(3) antiferromagnet for sufficiently large on-site repul- 
sions U. 

Motivated by both this prospect as well as the math- 
ematical challenges inherent to the problem, we propose 
several exact models for SU(3) spin chains in this Letter. 
The models are similar in spirit to the MG or the AKLT 
model for SU(2), and consist of parent Hamiltonians and 
their exact ground states. 

Consider a chain with N lattice sites, where N has to 
be divisible by three, and periodic boundary conditions 
(PBCs). On each lattice site we place an SU(3) spin 
which transforms under the fundamental representation 
3 = (1,0), i.e., the spin can take the values (or colors) 



blue (b), red (r), or green (g). (We label the represen- 
tations of SU(3) by their dimensions (the bold numbers) 
or their Dynkin coordinates (a pair of non-negative inte- 
gers) [33].) The trimer states are obtained by requiring 
the spins on each three neighboring sites to form an SU(3) 
singlet 1 = (0, 0), which we call a trimer and sketch it by 
o — o — o . The three linearly independent trimer states 
are given by 

and two more states obtained by shifting this one by one 
or two lattice sites, respectively. Introducing operators 
c\ a which create a fermion of color a (a = b, r, g) at 
lattice site i, the trimer states can be written as 



^trimer / 



. (a,/3, 7 )= 
(-^integer] 7r(b,r,g) 

(2) 

where fi = 1,2,3 labels the three degenerate ground 
states, i runs over the lattice sites subject to the con- 
straint that is integer, and the sum extends over all 
six permutations tt of the three colors b, r, and g. 

The SU(3) generators at each lattice site i are defined 

as 



x = \ E 



= 1,.... 



(3) 



cr,(j'— b,r,g 



where the A a are the Gell-Mann matrices [33]. The op- 
erators (3) satisfy the commutation relations [J?, J|] = 
5ij f abc Jf, a, b, c = 1, . . . , 8, with f abc the structure con- 
stants of SU(3). We further introduce the total SU(3) 
spin of v neighboring sites i, . . . , i + v — 1, 
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(4) 



where J \ is the eight-dimensional vector formed by its 
components (3). The parent Hamiltonian for the trimer 
states (2) is given by 



N , 
i=l v 
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To verify this Hamiltonian, note that since the spins on 
the individual sites transform under the fundamental rep- 
resentation 3, the SU(3) content of four sites is 

3 3® 3® 3 = 3-3 © 2-6 © 3-15 © 15', (6) 

i.e., we obtain representations 3, 6 = (0, 2), and the two 
non-equivalent representations 15 = (2, 1) and 15' = 
(4,0). All these representations can be distinguished by 
their eigenvalues of the quadratic Casimir operator [33] . 

For the trimer states (1), the situation simplifies as we 
only have the two possibilities 
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FIG. 1: Dispersions of the rep. 3 (left) and the rep. 3 trial 
states (right) in comparison to the exact excitation energies 
of (5). The lines are a guide to the eye. 



(where 3 = (0,1)), i.e., the total SU(3) spin on four 
neighboring sites can only transform under representa- 
tions 3 or 6. The eigenvalues of the quadratic Casimir 
operator for these representations are | and respec- 
tively. The auxiliary operators 



»-(W-|)((^)'-f) 



(7) 



hence annihilate the trimer states for all values of i, while 
they yield positive eigenvalues for 15 or 15', i.e., all other 
states. Summing Hi over all lattice sites i yields (5). 

There are two different kinds of domain walls between 
the degenerate ground states. The first kind consists of 
an individual SU(3) spin, which transforms under rep- 
resentation 3; the second kind consists of two antisym- 
metrically coupled spins on two neighboring sites, and 
hence transforms under representation 3. Since they can 
decay into each other, only one of these domain walls 
can constitute an approximate eigenstate of the trimer 
model. We have performed numerical studies on chains 
with N — 13 and N = 14, which clearly indicate that the 
elementary excitations of the trimer chain (5) transform 
under 3 (see Fig. 1). This result appears to be a general 
feature of rep. 3 SU(3) spin chains, as it was recently 
shown explicitly to hold for the HSM as well [16]. The 
elementary excitations of the trimer chain are deconfined, 
meaning that the energy of two localized representation 3 
domain walls or colorons does not depend on the distance 
between them. 

We now introduce a family of exactly soluble valence 
bond models for SU(3) chains of various spin representa- 
tions of the SU(3) spins at each lattice site. To formulate 
these models, we will use SU(3) Schwinger bosons b, 
(blue), r,r^ (red), and g,g^ (green) [35], which are de- 
fined by |b) = 4|0) = &t |o), |r) = c\ | } = r+ | ), 
and |g) = cj I } = 5 f | }, and satisfy [b, &t] = [ r , r+] = 
\_9i 9^] — 1 while all other commutators vanish. The 
Schwinger bosons can be used to combine spins trans- 
forming under the fundamental representation 3 = (1,0) 
symmetrically, and hence to construct representations 
like 6 = (2,0) and 10 = (3,0). 
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The trimer states (2) can be rewritten using SU(3) 

t r 1 „tl 



Schwinger bosons as IV'triLor) 



^^[6+,rT, <?Tj | ) with 



*"[6t,rt, fl t] 



JJ ( Y, signWal/?l +l7 ^ 2 



ti- u i ^ ( a '/ 3 -T)= 
(-3^ integer J 7r(6,r,g) 



(8) 

We obtain a representations 6 VBS from two trimer 
states by projecting the tensor product of two funda- 
mental representations 3 onto the symmetric subspace, 
i.e., onto the 6 in the tensor product 3 <£> 3 = 3 © 6. 
Graphically, we illustrate this as follows: 

one site 



T 

projection onto rep. 6 



(2,0) 



(9) 



This construction yields three linearly independent 6 
VBS states, which are readily written out using (8), 



(aO \ 

6 VBS / 



*"[6t, r t, fl t] 



|0) (10) 



for /i = 1, 2, or 3. These states are zero-energy ground 



states of the parent Hamiltonian Hsvbs = SjLi #» 

(ii) 

(see [34]). Note that the operators Jf , a = 1, . . . , 8, are 
now given by 6 x 6 matrices, as the Gell-Mann matrices 
only provide the generators (3) of the fundamental rep- 
resentation 3. As in the trimer model, two distinct types 
of domain walls exist, which transform according to reps. 
3 and 3. Since both excitations are merely domain walls 
between different ground states, there is no confinement 
between them. 

Let us now turn to the 10 VBS chain. By combining 
the three different trimer states (8) symmetrically, 

IVio vbsH* 1 [b\ r\ ff t] .$2 [ fe t r t 5 t] . tf» [ 5 t r t ff t] | o ) , 

(12) 

we automatically project out the rep. 10 in the tensor 
product 3 <g> 3 <S> 3 = 1 © 2 • 8 © 10 generated on each 
lattice site by the three trimer chains. This construction 
yields a unique state: 

one site 



(13) 




projection onto 10 = (3, 0) 



The parent Hamiltonian acts on pairs of neighboring 
sites. It is constructed by noting that the only rep- 
resentations which are included in both 10 © 10 and 
3 ® 3 © 3 © 3 (which is the representation content of the 
total spin of two neighboring sites of the VBS state as in- 
dicated by the dashed box above) are the reps. 10 = (0, 3) 



and 27 = (2, 2). With the eigenvalues of the Casimir op- 
erator, which are 6 and 8, respectively, we obtain the 
parent Hamiltonian for (12) 



N 



#io vbs = (( J * J m) 2 + 5 JiJi+i + 6) . (14) 

i=l 

The Hamiltonian (14) provides the equivalent of the 
AKLT model [6], whose unique ground state is con- 
structed from dimer states by projection onto spin 1, for 
SU(3) spin chains. 

Since the 10 VBS state (12) is unique, domain walls 
connecting different ground states do not exist. We hence 
expect the coloron and anti-coloron excitations to be con- 
fined in pairs, as illustrated below. The state between the 
excitations is no longer annihilated by (14), as there are 
pairs of neighboring sites containing representations dif- 
ferent from 10 and 27, as indicated by the dotted box 
below. As the number of such pairs increases linearly 
with the distance between the excitation, the confine- 
ment potential depends linearly on this distance, 
coloron anti-coloron 
3 3 




energy cost oc distance 

(15) 

The confinement force between the pair induces a linear 
oscillator potential for the relative motion of the con- 
stituents. The zero-point energy of this oscillator gives 
rise to a Haldane-type energy gap (see [36] for a similar 
discussion in the two- leg Hciscnbcrg ladder). We expect 
this gap to be a generic feature of rep. 10 spin chains 
with short-range antiferromagnetic interactions. 

Finally, we construct a representation 8 VBS state, 
where 8 = (1, 1) is the adjoint representation of SU(3). 
Consider first a chain with alternating reps. 3 and 3 on 
neighboring sites, which we combine into singlets. This 
can be done in two ways, yielding the two states 

o— O o—O o— O and -o O—o O-^J O . 
3 3 

We then combine one 3-3 state with the one shifted by 
one lattice spacing. This yields representations 3 © 3 = 
1©8 at each site. The 8 VBS state is obtained by project- 
ing onto the adjoint representations 8. Corresponding to 
the two 3-3 states illustrated above, we obtain two lin- 
early independent 8 VBS states, X F L and 5 ,R , which may 

be visualized as 
one site 



o—O 
-O o—O 



o—O 



o—O 
o- 



and 



-o O — o O — o O- 
O — o O — o O — o 



projection onto 8 = (1, 1) 



(16) 



These states transform into each other under parity or 
color conjugation (interchange of 3 and 3). The corre- 
sponding states may be formulated as a matrix prod- 
uct [34]. 
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The parent Hamiltonian for these states is constructed 
along the same lines as above, yielding 

#8vbs = E (V* J m) 2 + \ JiJi+i + I) ■ (17) 
»=i V - / 

^ L and \I/ R are the only (zero-energy) ground states of 
(17) for N > 3. 

The low-energy excitations of the 8 VBS model are 
given by coloron-anti-coloron bound states: 

anti-coloron coloron 
3 3 

o — O o— O c-<3 o—O o-O c^J 

c — O • O — o O — o • c — O 

* L * L (18) 

energy cost oc distance 

We find a linear confinement potential between the exci- 
tations, and hence a Haldane-type gap in the spectrum. 
Numerical studies on a chain with N — 8 sites provide 
evidence in support of this conclusion [34] . In addition to 
the bound state (18), the model allows for domain wall 
between the two ground states \E' L and SS? R . They con- 
sist of bound states of cither two colorons or two anti- 
colorons, which arc confined through the same mecha- 
nisms as the coloron-anti-coloron bound state (18), as 



one may easily infer from a cartoon similar to the one 
above. We hence expect a Haldane gap for each individ- 
ual domain wall as well. 

The results regarding confinement and dcconfincmcnt 
of the excitations of the 6, 10, and 8 VBS presented here 
are consistent with a rigorous theorem by Affleck and 
Lieb [37] on the existence of energy gaps in the spectrum 
of SU(n) nearest-neighbor Heisenberg spin chains. 

In conclusion, we have formulated several exact models 
of SU(3) spin chains. We first introduced a trimer model 
and presented evidence that the elementary excitations 
of the model transform under the SU(3) representations 
conjugate to the representation of the original spin on the 
chain. We further introduced three SU(3) valence bond 
solid chains with representation 6, 10, and 8, respec- 
tively, on each lattice site. The elementary excitations of 
the 10 and the 8 valence bond solid chain were found to 
be confined. 
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